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Abstract 

We develop a variational theory of hot nuclear matter in neutron stars and supernovae. It can 
also be used to study charged, hot nuclear matter which may be produced in heavy-ion collisions. 
This theory is a generalization of the variational theory of cold nuclear and neutron star matter 
based on realistic models of nuclear forces and pair correlation operators. The present approach 
uses microcanonical ensembles and the variational principle obeyed by the free energy. In this paper 
we show that the correlated states of the microcanonical ensemble at a given temperature T and 
density p can be orthonormalized preserving their diagonal matrix elements of the Hamiltonian. This 
allows for the minimization of the free energy without corrections from the nonorthogonality of the 
correlated basis states, similar to that of the ground state energy. Samples of the microcanonical 
ensemble can be used to study the response, and the neutrino luminosities and opacities of hot 
matter. We present methods to orthonormalize the correlated states that contribute to the response 
of hot matter. 

PACS numbers: 21.65. +f, 26.50. +c, 26.50.+X, 97.60. Jd, 97.60. Bw, 05.30.-d 
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I. INTRODUCTION 



Ah initio theories of strongly interacting hot matter are extremely challenging. In principle 
the properties of hot matter can be calculated starting from a realistic Hamiltonian with the 
path integral Monte Carlo method Calculations are practical for simple systems of 
interacting spin zero bosons such as atomic 4 He liquids and solids □ □ □ . They become 
more difficult even for simple systems of fermions interacting by spin independent potentials, 
such as atomic liquid 3 He |J, and hydrogen plasma jf| due to the fermion sign problem. The 
path integral Monte Carlo treatment is expected to become much more difficult due to the 
strong spin-isospin dependence of nuclear forces and their tensor and spin-orbit components. 
In the traditional Monte Carlo approaches these complexities of the nuclear forces make 
computations more expensive by a factor > 2 A , where A is the number of nucleons and 
the equality applies for pure neutron matter. With the present state of the art computing 
facilities traditional quantum Monte Carlo calculations have been carried out for cold neutron 
matter using a periodic box containing 14 neutrons jjj]. Attempts are also being made to 
eliminate this 2 A factor using the auxiliary field diffusion Monte Carlo method |sj], however 
the fermion sign problem is more acute for this method, and applications have been limited 
to cold pure neutron matter 

Cold nuclear matter has traditionally been studied with variational methods llCJ, lu| 

nn 

and Brueckner theory There is close agreement between these two methods, and 

comparison with the essentially exact Green's function Monte Carlo calculations suggests 
that the errors in present variational calculations of pure neutron matter are only ~ 8 % at 
densities < po = 0.16 fm -3 21] • In the case of symmetric nuclear matter the errors have been 
estimated to he < 10 % \14\. 

In this paper we develop the formalism for a variational theory for nuclear matter at 
finite temperature using correlated basis states (CBS) defined in the next subsection. The 
correlated basis states and the thermodynamic variational principle used to calculate the free 
energy of matter is discussed in the following subsections. In these subsections we review the 
scheme suggested in Ref. [Ha] to develop a variational theory of hot matter, and comment 
on the concerns expressed in its early applications till due to the nonorthogonality of 
the CBS. 
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In Section II we show that these problems can be resolved if one works in a microcanonical 
ensemble. We show that there are no orthogonality corrections to the free energy in this 
scheme. In Section III we consider the CBS that contribute to the response of the hot matter, 
and conclude in Section IV. 



A. Correlated Basis States 

Let the stationary states of a non interacting Fermi gas be denoted by |$j{nj(k, cr z )}), 
where {rij(k,a z )} are the occupation numbers of single particle states labeled with momen- 
tum k and spin projection a z , in the many-body state /. The single-nucleon states of a non 
interacting nucleon gas have isospin r z as an additional quantum number. We have sup- 
pressed it here for brevity. For each of the states /, we can construct a normalized correlated 
baS1S St ate (CBS) QQQ wMeh is e„— ly defined as: 

r.) , g|j> , J> a) 

VWWWT) 

where Q is a many-body correlation operator. Many problems in the variational theory of 
strongly interacting quantum liquids originate from the fact that useful forms of Q are not 
unitary operators. In recent studies Q has been approximated by a symmetrized product of 
pair correlation operators J 7 ^- [ill, llJ, |2I| where i and j label the nucleons: 

(2) 

i<j 

Here S stands for the symmetrization of the product of the pair correlation operators. In 
the present work we will assume this form of Q, however, improvements such as the inclusion 
of three-body correlations can be easily accommodated. The CBS obtained with this Q are 
not orthogonal to each other. The bras and kets with rounded parenthesis, ( | and | ) are 
used to denote these non orthogonal states, while the standard (| and |) imply orthonormal 
states. 

At zero temperature the parameters of Q or J-^ are determined variationally by minimizing 
the expectation value of the Hamiltonian, if, containing realistic interactions, in the ground 
state of the correlated basis |\I/o) obtained from the Fermi-gas ground state |$o)- The CBS 
are assumed to provide a good approximation for the stationary states of the interacting 
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system. Note that this is in accordance with an important assumption of the Landau theory 
of Fermi liquids, i.e. the stationary states (at least the low lying ones) of an interacting, 
normal Fermi liquid can be written in one to one correspondence with those of the non 
interacting one. 

If the no(k, a z ) be the occupation numbers of the single particle states in the ground state 
of the free Fermi gas, The |nj(k, a z ) — no(k, a z )\ can be interpreted as the quasi-particle 
{k > hp) and quasi-hole (k < k F ) occupation numbers of the CBS When the number 

of quasi-particles is finite the energy of the state /, Ej can be expressed as the sum of the 
ground state energy E , and a sum of quasi-particle and hole energies: 

Ej = =E + J^Mk, a z ) - n (k, ^)]e (k, a z ) . (3) 

k,cr z 

We assume that both the number of particles N and the volume of the liquid Q, go to oo at 
a fixed finite density p = N/fl. The density of quasi-particles and holes goes to zero when 
their number is finite. The single particle energies eo(k, a z ), have significant dependence on 
k near k F at low temperature, in addition to that absorbed in the 7^ term with m* 
being the effective mass of the quasiparticle. They are difficult to calculate ab initio. 

A correlated basis perturbation theory (CBPT) can be developed using the non orthogonal 
CBS Q, Q to study various properties of quantum liquids |2l|, [22| at zero temperature. 
Much later in the development of CBPT, a scheme to orthonormalize the CBS preserving 
their one to one correspondence with the Fermi gas states and the validity of Eq. (JBJ was 
found It simplifies CBPT considerably. 

The difference between the internal energies of a liquid at T > and T = is extensive, 
i.e. proportional to N, and thus infinite in the thermodynamic limit. This implies that at 
T > there is an extensive number of quasi-particle excitations, and the orthonormalization 
scheme of Ref. [3| can not be used without modifications. The present work can also 
be considered as an extension of that orthonormalization scheme to hot matter. At very 
small temperatures, the density of quasi-particles is small, and the T = formalism can be 
used neglecting the interaction between quasi-particles, as in Landau's theory. However, the 
domain of the applicability of that approach is very small 
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B. The Thermodynamic Variational Principle 



Let F(T) be the free energy of a quantum many body system at temperature T. All other 
arguments such as the density p and spin-isospin polarizations etc. have been suppressed for 



brevity. The Gibbs-Bogoliubov thermodynamic inequality |2J| states that 



F(T) < Tr (p v H) — TS V (T) , (4) 

where pv is any arbitrary density matrix (not to be confused with the density of the system 
p = N/fl) satisfying 

Trpy = 1 (5) 

and Sv(T) is the entropy of the density matrix py at temperature T. The equality holds 
when py is the true density matrix of the system. Typically pv is chosen to have the canonical 
form, 

exp(-(3H v ) 
Tvexp(-(jH v ) 

where (3 is the inverse temperature and Hy is chosen as a suitable, simple and variable 
variational Hamiltonian. In this case Eq. (J3J) becomes 

Tv(exp(- (3 H V )H) 
F[T) ~ Tr(exp(-/^)) " TSv{T) ' (7) 



The minimum value of 



TS V (T) , (8) 



Tv(exp(-(3H V )) 

obtained by varying Hy, provides an upper bound to the free-energy F(T). 

Schmidt and Pandharipande (Ref. [15], henceforth denoted by SP) proposed to use this 
variational principle to calculate properties of hot quantum liquids. They essentially ignored 
the nonorthogonality of the CBS and assumed that they are the eigenstates of Hy: 



#v|#i{ni(k,<r 2 )}) 



\^i{ ni (k,a z )}) . (9) 



The eigenvalues of this Hy can be varied by changing the single-particle energy spectrum 
€y(k, a z ) and the eigenfunctions by varying the correlation operator Q, or the pair correlation 
operators T%j. Note that the single particle energies depend on r 2 , p, T etc. but these 
dependencies are suppressed here. 



Hy has the spectrum of a one body Hamiltonian, since its eigenvalues depend only on 
the occupation numbers ni(k,a z ). It can therefore be easily solved. At temperature T the 
average occupation number of a single-particle state is given by 

n(k, o z = ±1) = emk J_^ ±) - i , (10) 

where the chemical potential \i± is required to satisfy 

/d k /_ v , , 

— n(k,* z = ±l). (11) 

In the above equation p± is the density of particles with a z = ±1. The entropy Sy(p,T) 
is given by j^J, 

S v (p } T) = -k B nJ2 J J^y N k '^) ln(n(k,a 2 )) + (l-n(k,(T,)) ln(l-n(k, a z ))] . (12) 

where ks is Boltzmann's constant. 

Since the CBS are not mutually orthogonal, Eq. (fT2|) is only an approximation if the 
variational Hamiltonian, H v , is defined by Eq. (JHJ). Eq. (JT2l) will be exact if orthonormalized 
correlated basis states (OCBS) are used instead of the non orthogonal CBS. If all the CBS are 

n 

orthonormalized by a democratic procedure (like the Lowdin method) |2||, which treats all 
the CBS equally, the diagonal matrix elements of the Hamiltonian H change by an extensive 
(oc N) quantity. 

The diagonal matrix elements of the Hamiltonian, H, can be evaluated using the standard 
techniques of cluster expansion and chain summation; these techniques have been developed 
and studied extensively in the variational theories of cold (zero temperature) quantum liq- 
uids. On the other hand, if all the CBS are orthonormalized using the democratic procedure 
mentioned above, the diagonal matrix elements of the Hamiltonian, H, in the corresponding 
OCBS are more difficult to evaluate systematically because of the extensive (oc N) orthog- 
onality corrections. As such, the variational theory of hot (finite temperature) quantum 
liquids so developed using the orthonormalization scheme discussed above, loses much of the 
simplicity of the corresponding zero temperature theory. 

At zero temperature a similar problem was addressed by identifying the ground state and 
the excitations about the ground state with a finite number of quasiparticles and quasiholes, 
as the important states which contribute to the equilibrium properties and linear response of 
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cold quantum liquids. It was shown in Ref. |2j| that a combination of democratic (Lowdin) 
and sequential (Gram-Schmidt) orthonormalization methods can be used to orthonormalize 
the CBS, such that the diagonal matrix elements of the Hamiltonian, if, are left unchanged, 
in the ground state and in the quasiparticle-quasihole excitations from the ground state. 

At a finite temperature, the many body states which contribute to the equilibrium prop- 
erties (free energy, specific heat etc.) and linear response of a quantum liquid are the 
many body states in the microcanonical ensemble at the corresponding temperature and 
the quasiparticle-quasihole excitations from them. (Zero temperature is a special case when 
the microcanonical ensemble consists of just one state viz. the ground state.) 

In this paper we will show that for a given (finite) temperature a statistically consistent 
microcanonical ensemble can be defined, such that when the CBS are orthonormalized us- 
ing a combination of democratic and sequential orthonormalization methods, the diagonal 
matrix elements of the Hamiltonian, if, are left unchanged for the many body states in the 
microcanonical ensemble and the quasiparticle-quasihole excitations from the microcanonical 
ensemble. This means that these matrix elements can be evaluated by borrowing methods 
directly from the zero temperature theory. 

As mentioned earlier, this work can be considered to be an extension of the orthonor- 
—„ sehetae of of Re, Q to finite tetaperatutes. However, it serves a m ore genera, 
purpose of introducing a variational theory at finite temperatures which has the same sim- 
plicity of formulation and efficiency in calculation as the corresponding zero temperature 
theory. 

II. VARIATIONAL THEORY IN A MICROCANONICAL ENSEMBLE 

In the previous section we have defined the non interacting Fermi gas states |$/) and 
the CBS Let us call the corresponding OCBS Note that the actual definition 

of l^j) will depend on how we choose to orthonormalize the CBS. We will denote any of 
these OCBS by and the actual orthogonalization procedure used to obtain them will 
hopefully be obvious from the context. Let us also define a 'microcanonical' subset .M(T), 
from the set of all labels I of the many body states (CBS, OCBS or non interacting Fermi 
gas) previously defined. We will call this set the 'microcanonical ensemble' at temperature 
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T. Henceforth the argument T will be supressed for brevity. Note that as yet we have not 
really said anything about which elements are included. We tackle this slightly non trivial 
problem in detail later in this section. For now, we assume that M. is a suitably defined 
'microcanonical' ensemble at the given temperature. We can legitimately define a density 
matrix, 



r ( 13 ) 



PMC =- . r 

NM ieM 

where Mm is the number of elements in the set Mi . 

It is well known in statistical mechanics that the thermodynamic averages of the den- 
sities of extensive quantities are the same in all ensembles: grand canonical, canonical or 
. In Eq. (JBJl we have used the canonical ensemble for the average value 

of H. 

With the microcanonical ensemble we obtain a simpler expression, 

<#> = T7- £ <*'l# l*'> • ( 14 ) 
Nm iTm 

In order to develop the variational theory of hot quantum liquids we have to orthonor- 
malize the CBS in the microcanonical ensemble, Ml. This can be easily achieved with the 
Lowdin transformation j^fj, 



2 

J e M 

+ l £ + ■ (15) 

J,K e M 

The coefficients 1, — |, |, • • • that occur in the Lowdin transformation are those which are 
found in the expansion of (1 + x)~ 1 / 2 . The overhead bar signifies, 



= ^>j\^ K ){l-5j K ) . (16) 

The orthonormal states are in one to one correspondence with the CBS and the 
Fermi-gas states |$r), and we are then justified in defining a variational Hamiltonian Hy 



such that 



#y|*j{nj(k,<7,)}> 



y^nj(k, a z )e v (k, a z 



.k,<r z 



l^jimik,^)}) , (17) 



EY\^j{ ni (k,a z )}) (18) 



thus removing the approximation inherent in Eq. The CBS G" M. are not orthonormalized 
by the transformation (fT5|) . Most of these states have little effect on the thermodynamic 
properties of the liquid in equilibrium at temperature T and density p. Formally these 
states should be first orthonormalized to those G M. by Gram-Schmidt's method, and then 
orthonormalized with each other using combinations of Gram-Schmidt and Lowdin methods 
This way their orthonormalization will have no effect on the states G M.. In the next 
section we will have occasion to discuss the orthogonalization of a subset of these states, 
viz. states with one (quasi)particle and one (quasi)hole with respect to the states in the 
microcanonical ensemble. 

In the variational estimate of the free energy [Eq. (J3J)] we should use the OCBS rather 
than the CBS. In the remaining part of this section we show that 

^i\h^j) = ^m^j) = , (19) 

i.e. if we define 



<5£j = - j)] (20) 

and 

5ej = ^5Ej (21) 

then, 

bei = (22) 

for / G M., in the limit iV — > oo. Therefore the variational free energy calculated with the 
SP scheme does not have any orthogonality corrections. 

At this point it is necessary to define the microcanonical ensemble (JM) more carefully. 
Typically, a microcanonical ensemble is defined as 

M = (J All states I with E MC < Ej < E MC + SE MC . (23) 

The actual value of 5E M q is unimportant as long as SE M q <C E M q. In our case it proves 
necessary that it takes a nonzero value. In the next subsection we will show that the simplest 
definition of A4, i.e. with 5E MC = 0, gives a divergent expression for the diagonal matrix ele- 
ments of H. This exercise will nevertheless help to illustrate some of the simplest elements of 
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the calculations that follow and will serve as a motivation for the following subsection where 
we formulate the problem slightly differently, which makes calculations more convenient, but 
is similar to defining Ai with a nonzero 5E M q. 



A. Energy Conserving Microcanonical Ensemble 

Let us define a set A4 , which we will call the Energy Conserving Microcanonical Ensemble 
(ECMC), as the set of all states with 

Ej = E MC (24) 

Consider a many body state / G A4 with, 

Ej = (tfj|if . (25) 

Then the change in the diagonal matrix elements of the Hamiltonian due to Lowdin or- 
thonormalization is given by 

= -~ J2 [^i\H-Ej\^j) (VjIV^ + ^jIH-E^j) +••• , (26) 

J e Mo 

where the dots denote higher order terms which can easily be obtained from Eq. (fTBTl . The 



nondiagonal CBS matrix elements (^j\H — Ej\^j) and (^//|^/j) can be evaluated with clus- 
ter expansions The leading two-body clusters contribute to the nondiagonal matrix 
elements only when two quasi-particles in J are different from those in I. Let quasiparticle 
states with momenta and k 2 be occupied in / and unoccupied in J, while states kj/ and 
k 2 ' be occupied in J but not in /. 

We will denote the CBS |*j) and the OCBS by 

= |k r ,k 2 , : J-k^ka) , (27) 
|*j) = |k r ,k a ,: J-k^ka) . (28) 

Note that in this notation can be written as |ki, k 2 : / — ki, k 2 ) 
The J — > J transition occurs via the scattering of two quasi-particles from states 
(ki,k 2 ) — > (ki/,k 2 /). Momentum conservation implies 

k x + k 2 = ki, + k 2 , . (29) 
10 



Unless this condition is satisfied, the nondiagonal CBS matrix elements are zero. 



The two-body cluster contributions to (^j\H — Ej\^fj) and are respectively given 

by 

(k^k^-k^k^^,^) and (k r ,k 2 , -k 2 ,,ki/|(^g -l)|ki,k 2 ) , (30) 
where the two-body effective interaction is given by: 

h 2 



% = Fa 



VijFi 3 - — (V%- + 2V^ y • V) 



(31) 



m 

the bare two-body interaction is denoted by v^, and the non interacting two-particle states 
are 

(r i ,r j |k 1 ,k 2 ) = ie J ( kl ^ +k ^) . (32) 

The factor 1/Q comes from the normalization of the plane waves. We have suppressed the 
spin wave functions for brevity. 

The CBS matrix elements can be represented by diagrams, such as those in Figs, d [2] and 
El which have been analyzed in detail in j^]. We will adopt their notation and use their 
results. In all the diagrams we use the following conventions. 

• The points in these diagrams denote positions of the particles: r i; r,-, .... 

• The dashed lines connecting points i and j represent correlations, i.e. terms originating 
from Ffj — 1. When T, L j = /(%) this notation is sufficient, however, a more elaborate 
notation for the correlation lines is needed when T is an operator with many terms 
j^J. For brevity we will show diagrams assuming Tij = /(r^), commonly called the 
Jastrow correlation function. 

• The solid lines represent v^J . There can only be one solid line in a diagram representing 
matrix elements of H. 

• The lines with one or two arrowheads represent state lines. The arrowheads are labeled 
with quasi-particle states. A state line with a single arrowhead labeled kg going from 
point i to point j indicates that the particle i is in state k^ in the ket and particle 
j is in k^ in the bra Diagrams representing diagonal CBS matrix elements can 
have state lines with only one arrowhead, since the state k^ is occupied (or unoccupied) 
in both the bra and the ket. 
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Diagrams contributing to the nondiagonal CBS matrix elements have state lines with 
two arrowheads. The number of these lines equals the number of quasi-particle states 
that are different in and A state line with arrowheads and k^, going 

from i to j indicates that i is in state k^ in the ket, while j is in state ki> in the bra, 
and that k^ and k^ are unoccupied in the bra and the ket respectively. 

Only one state line must emerge from a point and only one must end in a point because 
each particle occupies only one quasi-particle state in the bra and the ket. This implies 
that the state lines form continuous loops. In direct diagrams the state lines emerge 
from and end on the same particle, while in exchange diagrams they connect pairs of 
particles. 

• The contribution of a diagram is given by an integral over all the particle coordinates 
rj in the diagram. The integrand contains factors of (/ 2 (^) — 1) for each correlation 
line, for the interaction line, e lkrFl / \/VL for each state line t emerging from a point 
rj and e~ lki '' rj / VTi for each state line £' ending in rj. 

The two-body (26) direct (d) and exchange diagrams (e) representing the nondiagonal 
matrix elements (|30|) are shown in Fig. [TJ The contributions of the 2b.v e ^ diagrams are 
given by 



2b.v eff .d = ~J rf%exp(-^(k 1 , - k 2 ,) • ri > e// (r) exp^^ - k 2 ) • r -) , (33) 
2b.v eff .e = ~ J d 3 nj exp(-^(k 2 , - k v ) ■ r l3 )v eff (r) exp^k, - k 2 ) ■ r y ) , (34) 

for spin independent v (and /) . The contributions of 2b. T 2 diagrams are obtained by 
replacing v e ^ by T 2 — 1. All 2b diagrams have a contribution of order 1/fl. 

The change in energy, 5E l7 [Eq. (|26|) ] contains products of the 2b.v e ^ and 2b. T 2 diagrams. 
These products are of order 1/fl 2 . The total contribution of the leading 2b cluster terms 
to S Ei is obtained by summing over the states ki, k 2 , ki/, k 2 '. Each allowed combination of 
these states corresponds to a many-body state in the set M.q. The quasi-particle states k x 
and k 2 can be any two of those occupied in I. Thus the sum over these gives a factor of 
order N 2 . Next we sum over kj/ and k 2 /. The total momentum k^ + k 2 / is determined from 
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Eq. The magnitude of the relative momentum, 



(35) 



is constrained by energy conservation, 



ey(ki) + e y (k 2 ) = e v (k r ) + e v (k 2 >) . 



(36) 



required for J to be in the set M.q. Thus the sum over states allowed for k]/ and k 2 ' 
corresponds to an integration over the direction of k^/. It gives a factor of order fi 1 / 3 . 
Hence 



i.e. 5ej.2b — > as iV — > oo. Note that with the constraint of momentum conservation 
alone we can integrate over the magnitude of ki/ 2 '. This integration gives a factor of order 
fl and makes 8ej.2b of order 1. The equal energy constraint in the ECMC Aio ( Eq. J25)) 
makes 6ei.2b vanish in the thermodynamic limit. 

The above analysis can be carried out for contribution of clusters with three or more 
particles to 5Ej. Consider, for example, states J which differ from / in occupation numbers 
of three quasi particles. These states can be reached by scattering three quasi-particles in 
/ in states ki,k 2 ,k 3 to states ki/,k 2 /,k 3 / occupied in J. The relevant, direct 36 diagrams 
are shown in Fig. El Each diagram is of order 1/Q 2 . The contribution of three body cluster 
terms to 5Ej are products of v e ** and T 2 — 1 diagrams. Each of these products is of order 
1/Q A . We get a factor of order N 3 by summing over states k 1; k 2 and k 3 , and a factor 
f2 4 / 3 by summing over k]/,k2',k3/ with constraints of momentum and energy conservation. 
Thus their total contribution to Sei is of order A^^^fi -4 ^ 3 ^ 4 / 3 ) ~ p 2 fi -2 / 3 which vanishes 
in the thermodynamic limit just like the contribution from the leading two body cluster 
terms. Similarly the contribution from all connected terms can be shown to give vanishing 
contribution to bej in the thermodynamic limit. 

The terms of Eq. (|26|) will also contain disconnected diagrams like those shown in Fig. El 
These diagrams by themselves will give rise to unphysical divergent (non extensive) contri- 



5e T .2b 




(37) 



(38) 
(39) 
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butions to the energy. To extract any physically meaningful result from the theory, these 
diagrams must cancel identically. 

Disconnected diagrams in the expansion for the shift in energy (SEj) can be classified 
into two types. 

• Diagrams in which each connected cluster conserves energy. 

• Diagrams in which only the whole diagram conserves energy, i.e. each connected cluster 
does not conserve energy. 

Let us, consider the first case. 

Consider the simplest possible divergent diagrams, i.e. Fig. El 
Let us denote the CBS \9 K ) and \ty L ) by 

= |ki/,k 2 /,k 3 /,k4/ : I - k 1 ,k 2 ,k 3 ,k 4 ) 
\V K ) = Ikx^ky : J-ki.ka) 

|* £ ) = |k 3 ,,k 4 , : J-k 3 ,k 4 ) (40) 

Let ki, k 2 and ky, k 2 / conserve momentum and energy amongst themselves and similarly 
for k 3 , k 4 and k 3 /, k 4 /. 



ki + k 2 = k X / + k 2 / 
k 3 + k 4 = k 3 / + k 4 / 



(41) 
(42) 



and 



e y (k 1 ) + e y (k 2 ) = ev(k v ) + e y (k 2 ,) , (43) 
e y (k 3 ) + ey(k 4 ) = e y (k 3 ,) + e y (k 4 ,) . (44) 

An inspection of the series on the right hand side of Eq. ([2~K|) will show that thirteen (13) 
terms in total will give rise to an integral which will be represented by products of diagrams 
shown in Fig. El In Table |U we list the terms and also their corresponding prefactor in 
the series. As one can see, the sum of the prefactors is identically zero. Thus products 
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of diagrams of the type shown in Fig. El has no contribution to the change in energy per 
particle (Sei). The cancellation of corresponding exchange diagrams and all other divergent 
diagrams of this order with three or more body connected pieces can also be shown to cancel 
with analogous book keeping. The divergent diagrams of the next highest order can also be 
shown to cancel identically. 

Now consider the case when the individual clusters do not conserve energy, i.e. 
Eqs. (|1U H2J) are still true but Eqs. ((331 133} are not true. Instead, 

e v (ki) + e y (k 2 ) + ey(k 3 ) + e y (k 4 ) = e v (k v ) + e v (k 2 /) + e v (k 3 ,) + e v (k 4 ,) . (45) 

In this case the states K and L no longer belong to the same ECMC as I and J. Thus, 
none of the terms in Table 1, except for the first, will be included in the sum, i.e. the 
divergent terms, (^i\v e ^ fifj and its complex conjugate will not get cancelled. 
The total contribution of terms like these to SEj is of order p 4 f2 4//3 , i.e. the shift in the 
energy per particle, 5ej ~ p 3 f2 1//3 , diverges in the thermodynamic limit. 

The survival of divergent terms is rather artificial and arises from the sharp energy con- 
servation constraint that we imposed on the states in M. Q . This provides the motivation 
to define a ensemble where this constraint is relaxed slightly. In what follows we will show 
that this can be done consistently where none of the divergent terms are present while the 
diagonal matrix elements of the Hamiltonian are preserved. 



B. The set of Most Probable Distributions (MPD) 

Consider an ideal gas of fermions in a box of volume L 3 (= Q) . The single particle energy 
levels are given by 

h 2 n 2 

where m is the mass of the fermions and rij is a vector with integer components. In what 
follows we will use units where = 1. Let the total number of particles in the box be N. 
The density of states at any single particle energy ey is given by 

g(e v ) - L 3 ef (47) 

Now consider a cell with an energy width of around an energy level ey . Then the 

15 



number of single particle energy levels in this cell is 

u{e v )^g{e v )—r^Lh 1 ^^. (48) 



The exact value of A is not important, except for the fact that it is dimensionless and of 
order 1 or less. We can always choose A so that w(ey) is large, 

u{e v ) > 1 . (49) 

Let the total number of particles in the energy cell around ey be n{ey). Let S{n{ey)) be 
the entropy corresponding to the configuration (distribution) n(ey). It can be easily shown 
that the entropy, Sy corresponding to the most probable distribution of number of particles 
per energy cell , subject to the constraints 

5>(ey) = N (50) 
and ^e v n(e v ) = E MC , (51) 

is given by Eq. (fT2|) . i.e., the distribution n(ey) in addition to satisfying Eqs. (|50|51|) also 
satisfies the maximization condition, 

S(n{e v )) = Maximum(S(n(e y ))) (52) 

The temperature and the chemical potential are the Lagrange multipliers of the minimization 
procedure. Now we will define our microcanonical ensemble as the set A4 of all configurations 
whose cell distribution is n(ey). We will call this the set of Most Probable Distributions 
(MPD). Please note that this is different from defining an ensemble with total energy E MC , 
but is roughly the same as defining an ensemble with an average energy E MC and a small 
non zero energy width 5E MC . 

M. = ^J( All states with the distribution n(ey)) (53) 

It should be emphasized here that none of the conclusions that follow depend explicitly on 
the actual single particle spectrum given by Eq. (j4*H| . All the conclusions remain unchanged 
as long as the single particle energy levels produce a continuum in the thermodynamic limit. 
However we will continue to use Eq. (|46|) because the calculations are more transparent this 
way. 
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The probability of fluctuations about the most probable distribution is given by Einstein's 
relation, 

P(5n) ~ e~ ss , (54) 

where P(Sn) is the probability of a fluctuation of size 5n and 5S is the corresponding decrease 
in entropy. Around the most probable distribution, 

5S ~ 5n 2 , (55) 

i.e. the probability of fluctuations vanishes exponentially with the size of the fluctuations. 

In addition, the fluctuation (standard deviation) in the value of the total energy, Ef , in 
A4, can be easily shown to be, 

5E < (56) 

i.e. 5E is non macroscopic; the fluctuation in the energy per particle vanishes in the ther- 
modynamic limit. Thus, the ensemble we have defined is a consistent one in the statistical 
sense. 

Now let us discuss the allowed scattering processes within the set M.. For two states to 
be in A4, they must have the same populations n(ey) in all the energy bins (cells). What 
this means is that they must be connected to each other through excitations within cells. 
For example let / and k]/, k 2 / : / — ki, k 2 , be elements of M.. Let 

ev(ki) < e y (k 2 ) 
and ey(ki') < ey(k 2 ') 

Then, we need to have 

|ey(ky) - Cy(ki)| < 
and |ey(k 2 ') - ey(k 2 )| < 

where and are the widths of the cells containing k x and k 2 respectively. Please note 
that this, in general, means that there is no exact energy conservation, but that there is 
approximate energy conservation for each individual particle. 



Ai 
L 2 
A2 

L 2 



(57) 
(58) 



(59) 
(60) 
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Let us discuss the orthogonalization correction in M.. Consider the disconnected diagrams 
first. Following the discussion before, let us define states J, K and L as in Eq. ijlHJ) . with, 

ey(ki) < e v (k 2 ) < e v (k 3 ) < e v (k 4 ) , 
e v {k v ) < e v (k a ,) < e v (k 3 ,) < e v (k 4 ,) . (61) 

Let us assume that they obey Eqs. pi| H2|) i.e. they form two disconnected momentum 
conserving clusters. 

Now it is crucial to observe that if / and J belong to M. then ki and ky must belong 
to the same energy cell; similarly for k 2 and k 2 ', k 3 and k 3 / , and k 4 and k 4 /. This implies 
that K and L must also belong to M.. This was not the case when we had merely imposed 
overall exact energy conservation. 

Thus, in this case all the terms in Table |U will contribute to the sum in Eq. (|26|) . As such 
the disconnected diagrams will cancel each other and we will be left with a connected, non 
divergent sum, i.e. 

^disconnected = q ^ 

identically. The problem with disconnected diagrams that we encounter in ECMC is resolved 
in MPD. 

For the sake of completeness, we show that the connected diagrams also have a vanishing 
contribution towards the energy in the thermodynamic limit. Consider the two body cluster 
contributions to 5Ej. The 5Ej [Eq. (|26|) ] contains products of the 2b.v e ^ and 2b. T 2 diagrams. 
These products are of order 1/Q 2 . The total contribution of the leading 2b cluster terms to 
SEi is obtained by summing over the states ki, k 2 , ki' and k 2 '. Each allowed combination 
of these states corresponds to a many-body state in the set M.. The quasi-particle states k x 
and k 2 can be any two of those occupied in /. Thus, the sum over k x and k 2 gives a factor 
of order N 2 . Next we sum over and k 2 /. The total momentum + k 2 ' is determined 
from Eq. (|29|) . The magnitude of the relative momentum: 

k V2 , = l(k; - k' 2 ) , (63) 

is constrained by Eqs. l6f)|) . 

The sum over states allowed for kr and k 2 > corresponds to an integration over ki/ 2 '. But 
ki' 2 / is constrained to lie in a shell of width A where, A ~ min(Ai, A 2 ). Thus, the sum over 
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ki' and k 2 ' gives a contribution ~ fi-p up to a factor of order 1 (the factor of f2 comes from 
the density of states). Therefore the total contribution of the 2b diagrams after summing 
over ki, k 2 , k^, k 2 ' is, 

5Ej.2b ~ Vlp 2 ^ . (64) 
Thus the shift in the energy per particle is, 

6ej.2b ~ , (65) 

which vanishes in the thermodynamic limit. 

The above analysis can be easily carried out for contribution of connected clusters with 
three or more particles to the 5Ej. Consider for example states J which differ from I 
in occupation numbers of three quasi particles. These states can be reached by scattering 
three quasi-particles in I in states k 1; k 2 , k 3 to states k x /, k 2 ', k 3 / occupied in J. For example 
consider the direct 3b term shown in Fig. EJ. Each is of order l/fi 2 , thus their contribution to 
5Ej is of order l/f2 4 . We get a factor of order N 3 by summing over k x , k 2 , k 3 , and a factor 
Q 2 (pr) 2 by summing over ki/,k 2 /,k 3 / with constraints of momentum and (approximate) 
energy conservation. Thus their total is of order Qp 3 (j^) 2 - The contribution to the shift in 
energy per particle is 5ej ~ p 2 (77) 2 which vanishes in the thermodynamic limit: similarly 
for higher order clusters. Therefore, for connected clusters we see that 

^connected _^ q ^ ^gg^ 

in the thermodynamic limit. 

Thus, as claimed earlier in the section we have shown that it is possible to define a 
statistically consistent microcanonical ensemble such that Eq. (|T9|) is true for the elements 
in the microcanonical ensemble Ai, 

C. Discussion 

The simplest choice for a microcanonical ensemble is the ECMC. The ECMC has the 
following properties: 

• The states in ECMC have exact energy conservation; this imposes a sharp energy 
cutoff. 
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• Arbitrarily high single particle energy transfers are allowed while still remaining in the 
same ECMC 

It is due to the second property that clusters in disconnected diagrams can have arbitrary 
energy transfers and hence divergent contributions. These contributions are normally (in 
a canonical ensemble, when all states are included) canceled by contributions from higher 
order terms. But by imposing exact energy conservation we exclude the states which lead to 
these higher order terms which cancel the divergent part. Thus, we are left with a divergent 
series. 

In MPD on the one hand we relax the energy conservation slightly, and on the other hand 
we limit single particle energy transfers to the width of the energy cells. We showed that 
this leads to a convergent series. Also, we showed that the total energy is well defined in 
a MPD and that states with large deviations from the MPD (i.e. states with large single 
particle energy transfers) are exponentially improbable. 

The main difference between ECMC and MPD is that one follows from a conservation 
law and the other from a distribution. There can be states in the ECMC whose population 
distribution in the energy cells is very different from the MPD n(ey), but as long as the total 
energy of the state, E v = E MC , this state is a valid member of the ECMC. However, the 
number of these states is negligible as compared to the total number of states which have 
the MPD, and hence they can be safely neglected in the thermodynamic limit. On the other 
hand, those states whose total energy Ey differ from Emc by a non macroscopic amount and 
have the same population distribution as the MPD should be included in the microcanonical 
ensemble. We have shown that, for our purposes, a typical state in a microcanonical ensemble 
is given by the (most probable) distribution and not by an exact conservation law. 

Thus, we have shown that a consistent choice for M. does exist. We have also shown that 
the most obvious choice, namely, the ECMC leads to divergences in the theory. We traced 
these divergences to the existence of the sharp cutoff due to the exact energy conservation 
imposed on the states. Then we showed that these unphysical divergences can be removed by 
relaxing the energy conservation slightly, with the set of MPDs. We showed that MPDs can 
be consistently treated as microcanonical ensembles and that the diagonal matrix elements 
of the Hamiltonian remain unchanged upon orthogonalization in this ensemble. 
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In practical calculations, a microcanonical sample, M. of CBS is given by 



\^{n MC (k, a z )}) 



(67) 
(68) 



riMc(k, <J Z ) = 1 with probability n(k, o~ z ); else zero . 




(69) 



Since the Hamiltonian Hy can be easily solved, we can find the temperature corresponding 
to this energy. In the N — > oo limit it is just that used to find the n(k, er 2 ) [Eq. (fT0|) ]. All the 
MC states belonging to this set can be found by allowing particles in \^mc) to scatter into 
allowed final states. Each scattering produces a new CBS belonging to the same MC set. We 
denote this set by M.. If the quantum liquid is contained in a thin container with negligible 
specific heat, then it passes through the states in M. when in equilibrium at temperature T 
and density p. 

Eqs. (ffiTl l68j) have been recently used to calculate the rates of weak interactions in hot 
nuclear matter Note that \^mc) is a CBS since n MC (k, u z ) are either 1 or 0. When 

the number of particles in \^mc) is large the fluctuations due to sampling the probability 
distribution n(k, a z ) are negligible, and this state has the desired densities p± and energy per 
particle appropriate for the desired temperature T and Hamiltonian Hy used to calculate 
the n. Neglecting these fluctuations in the limit iV — > oo we obtain the variational estimate 
for the free energy, 



where the minimum value is obtained by varying the Q and ey(k, a z ). The {^mc\H\^mc) 
can be calculated with standard cluster expansion and chain summation methods used in 
variational theories of cold quantum liquids Q- At low temperatures (<C T F ) this 
method is particularly simple because the zero temperature Q and e(k, a z ) provide ve ry g ood 
approximations to the optimum. The main concerns raised in past applications 
of the SP scheme is that it neglects the nonorthogonality of the CBS, and provides only 
upperbounds for the free energy. Here we address only the first. At zero temperature the 
difference between the variational and the exact ground state energy has been estimated 




(70) 
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with correlated basis perturbation theory [21]. It may be possible to extend these methods 
to finite temperatures. 



III. ORTHONORMALIZATION OF THE QUASIPARTICLE-QUASIHOLE EXCI- 
TATIONS 

In the calculation of nuclear response functions one needs to use the diagonal matrix 
elements of the Hamiltonian in the quasiparticle-quasihole states [22]. At least at zero 
temperature the leading contribution to the dynamic structure function comes the lp-lh 
states. Here we will limit our discussion to the diagonal matrix elements of the Hamiltonian 
in the lp-lh excitations from the states in M.. 

Consider a OCBS I E. JA, where the single quasiparticle state with momentum 

h is occupied but the single quasiparticle state with momentum h + k is not. We will 
denote the quasiparticle-quasihole OCBS where the quasiparticle state with momentum h 
is replaced by one with momentum h + kby |h + k: J — h), and the corrresponding CBS 
by |h + k : I — h). The CBS |h + k : / — h) is orthogonal to all the states in A4, because 
they have different total momenta. Thus it only needs to be orthonormalized with all the 
other quasiparticle-quasihole excitations with the same momentum, via the Lowdin method. 
The excitations with two or more quasiparticles and quasiholes should be orthonormalized 
with the quasiparticle-quasihole states using a sequential method. But this does not have 
any effect on the quasiparticle-quasihole states, so we do not discuss them any further. 

The quasiparticle-quasihole OCBS is given by 

|h + k : J-h) = |h + k : I-h)-- ^ |h' + k : I'-ti) (h' + k : I' - h'|h + k : / - h) + - • • , 

I'eM 

(71) 

where h' G I'. The diagonal matrix elements of H are given by 

(h + k : / - h|#|h + k : / - h) = (h + k : / - h|#|h + k : I - h) 

- - t( h + k : 1 - h \ H \ h ' + k : /' - h') 

I'eM 

x (h' + k:/'-h , |h + k:/-h) + c.c.] + --- (72) 
In actual calculation of response functions one needs the difference between (h + k : 



22 



I - h|#|h + k : / - h) and Let us define 

E ph = [(h + k : I — h|#|h + k : / - h) - (^1^1^)] . (73) 

At zero temperature, in accordance with Landau's theory, one can define single particle 
energies for quasiparticle and quasiholes as done in Eq. [BJ The quantity E p h is analogous to 
(i.e. is a finite temperature generalization of) the difference between the quasiparticle energy, 
eo(h + k), and the quasihole energy, eo(h). We will show that E p h has no orthogonality 
corrections. 

The orthogonality correction to E p h is given by 

5E ph = [(h+ k : / - h|#|h + k : / - h) - (tf jI^I*/)] 

- [(h + k : I - h|i?|h + k : I - h) - (^1^1*/)] (74) 

-J2 [(h + k : / - % e// |h' + k : /' - h') (h' + k : /' - h'|h + k : / - h) 
. I'eM 
+ c.c.l + - ■ ■ 1 



-J2 [( h:I ~ % e// |h' : /' - h') (h' : /' - h'|h : / - h) 

. I'eM 

+ c.c.] + ---] (75) 



The non diagonal matrix elements of the Hamiltonian (v e ^) and unity in the second equation 
are of order 1/Q or less. 

The shift SE ph will contain both connected and disconnected terms. The disconnected 
terms can be shown to cancel exactly using arguments similar to the ones used in the last 
section. We will consider the connected terms only. 

The state |h' + k : I' — h') can be of the following types: 

• Type 1 : h' = h , J ' ^ I 

• Type 2 : h' ^ h. 

For Type 1 terms, the terms of the matrix elements (h + k : / — h|t> e ^|h + k : /' — h') 
and (h + k : /' — h|h + k : / — h) do not depend on k. (The contribution of the terms which 
contain the exchange line h + k vanish in the thermodynamic limit as compared to the 
leading order terms. This can be easily seen by explicitly writing down the cluster expansion 
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for the matrix elements.) Thus the contribution of these matrix elements is canceled by the 
corresponding terms (h : I — h|t> e -^|h : V — h) and (h : I' — h|h : I — h). 

Consider the Type 2 CBS. Let hi be the single quasiparticle state in I' which is in the 
same energy cell (as defined in the last section) as h, and let h' x be the single quasiparticle 
state in / which is in the same energy cell as h'. Since both / and I' belong to A4, there 
will be at least one choice for hi and h' 1? although in general their choice is not unique. 

|h' + k : /' - h') = |h' + k, hi : /' - h', h 1 ) . (76) 

Similarly |h + k : / — h) can be written as |h + k, h' x : / — h, hj). For Type 2 CBS the 
leading contribution to Eq. (|74*Jl comes from states where V — h', hi = I — h, h' x . In 
this case each of the matrix elements (h + k, h'i : / — h, h'-Jv^lh' + k, hi : V — h', hi) and 
(h' + k, hi : /' — h', hi|h + k, h' x : I — h, h' x ) are 1/fl. Also, h^ can be any of the occupied 
quasiparticle states in I. Hence summing over h' x gives a factor of N. The sum over h' 
and hi along with momentum conservation and approximate energy conservation gives a 
term of order Note that there is no summation over h. Thus the total leading order 

contribution from the Type 2 states is of order p-p-, which vanishes in the thermodynamic 
limit. 
Thus, 

5 E ph -> . (77) 

There is no orthogonality correction to the energy differences which enter the calculations 
of response functions. 

IV. CONCLUSION 

We have developed a variational theory for hot quantum liquids. We have shown that 
the correlated basis states which provide a reasonable description of the ground state of 
quantum liquids can be used to describe quantum liquids at finite temperature. Although 
the correlated basis states are not orthogonal to each other by construction, the free energy 
calculated in a suitably defined microcanonical ensemble of the correlated basis does not 
have any corrections due to nonorthogonality. As such the powerful cluster expansion and 
chain summation methods developed for zero temperature quantum liquids can be used at 
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finite temperature without any reformulation. We have also shown that the energy differences 
which are needed for calculating response functions do not need any orthogonality corrections 
either. 

We wish to emphasize that the arguments used in this work do not depend on the detailed 
form of the correlation functions or the choice of the trial single quasiparticle spectrum. 
The correlated functions are merely required to be sufficiently well behaved so that all the 
integrals used in the Sections II and III are finite. Any reasonable form for the correlation 
functions can be expected to satisfy this requirement. Although we have used only two 
body correlation functions without any state dependence or backflow terms to illustrate our 
results, the arguments can be easily extended to include both of the above and also three 
body correlations. Similarly, the arguments can also be extended for any trial single particle 
spectrum which has a vanishing energy level spacing in the thermodynamic limit. 

We have not addressed the problem of non diagonal matrix elements, which are required 
for certain applications. One such case is presented in the calculation of weak interaction 
rates, where the relevant matrix elements are the non diagonal matrix elements of one body 
operators. Work is in progress to calculate the non diagonal matrix elements of these one 
body weak interaction operators including the orthogonality corrections. 

We have also not discussed the actual forms for the correlation functions or the trial 
single quasiparticle spectrum which may be useful for calculating the free energy at finite 
temperature. At low enough temperatures the zero temperature forms may provide good 
approximations, but at higher temperature this is probably not true. Methods to optimize 
free energy computations at finite temperature are being developed. 
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(tf/lt^ltfjKtfjl*/) 

(*j|« e ^|*A-)(*«:|*L)(*L|*/) 

(ttili^/ltfjJOPjI^a-Xtfjrltfj) 
(*j|t; e// |*ji:)(*jf|*j)(*j|*j) 





u e// |*j)(*j|*j) 
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u e// |*j)(*j|*z)(*i|*j) 
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3 
16 


(*/|* L )(*L 


*j)(*j|« e// |*jc)(*jf|*/) 


3 
16 


Total 








Table I: The contribution of all the terms which give rise to a term like Fig.0 
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Figure 1: All two body cluster diaj 
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Figure 2: Examples of three body connected diagrams 
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Figure 3: Examples of disconnected diagrams 
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